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ABSTRACT 

Gamma-ray bursts (GRBs) are proposed as candidate sources of ultra-high energy cosmic rays (UHECRs). 

We study the possibility that the PeV neutrinos recently observed by IceCube are produced by GRB cosmic 
rays interacting with the interstellar gas in the host galaxies. By studying the relation between the X-ray ab¬ 
sorption column density Ah and the surface star-formation rate of GRB host galaxies, we find that Ah is a good 
indicator of the surface gas density of the host galaxies. Then we are able to calculate the neutrino production 
efficiency of CRs for GRBs with known Ah- We collect a sample of GRBs that have both measurements of Ah 
and accurate gamma-ray fluence, and attempt to calculate the accumulated neutrino flux based on the current 
knowledge about GRBs and their host galaxies. When the CR intensity produced by GRBs is normalized with 
the observed UHECR flux above ~ lO'^eV, the accumulated neutrino flux at PeV energies is estimated to be 
about (0.3 ± 0.2) x 10“^GeV cm“^ s“* sr“' (per flavor) under the assumption that GRB energy production rate 
follows the cosmic star-formation rate and the favorable assumption about the CR diffusion coefficient. This 
flux is insufficient to account for the IceCube observations, but the estimate suffers from some assumptions in 
the calculation and thus we can not rule out this scenario at present. 

Subject headings: neutrinos - gamma-ray; bursts 

1. INTRODUCTION 

Gamma-ray bursts (GRBs) have been proposed as a poten¬ 
tial origin of ultra-high-energy cosmic rays (UHECRs) (e.g., 

Waxman 1995; Vietri 1995; Milgrom and Usov 1995). High- 
energy neutrinos are thought to be a useful messenger to probe 
CR acceleration in GRBs, as they are predicted to be pro¬ 
duced in the dissipative fireballs, where cosmic-ray protons 
interact with fireball gamma-ray photons through photopion 
process (e.g., Waxman & Bahcall 1997; Guetta et al. 2004; 

Dermer & Atoyan 2006; Murase et al. 2006; Wang& Dai 
2009). However, search for high-energy neutrinos in coin¬ 
cident with GRBs using IceCube has failed to find any as¬ 
sociated neutrinos so far (Abbasi et al. 2012; Aartsen et al. 

2015). The non-detection has put stringent constraints on the 
neutrino production efficiency and fireball properties of GRBs 
(He et al. 2012; Zhang & Kumar 2013; Gao et al. 2013; Ya- 
cobi et al. 2014). Because the low neutrino flux could be 
simply due to a low efficiency in converting CRs to neutrinos, 
as expected in some GRB models with large dissipation ra¬ 
dius, one can not rule out CR acceleration in GRBs with the 
non-detection of neutrinos. 

CRs accelerated by GRBs will finally escape from the 
source and enter the interstellar space of the host galaxy. 

The proton-proton collisions between CRs and nuclei in the 
interstellar medium will also produce high-energy neutri¬ 
nos. Recently, IceCube Collaboration reported the detec¬ 
tion of extraterrestrial TeV-PeV neutrinos with a best-fit flux 
of = (0.95 ± 0.3) X 10“^GeV cm“^ s“' sr“* per flavor 
(Aartsen et al. 2014). Noting the coincidence between 
the IceCube neutrino flux and the Waxman-Bahcall bound 
(^ 10“*GeVcm“^s“'sr“') derived from the flux of UHECRs 
above lO'^eV, some authors suggest that IceCube neutrinos 
may be produced by the same source responsible these UHE¬ 
CRs (Waxman 2013; Katz et al. 2013). As GRBs are one 
candidate source of UHECRs, one may wonder whether neu¬ 
trinos resulted from the CR collisions with the ISM in the host 


galaxy can explain the IceCube observations (Waxman 2013, 
Wang et al. 2014). In this paper, we study this interesting pos¬ 
sibility by perform a calculation of the expected neutrino flux 
with the knowledge about the properties of GRBs and their 
host galaxies. 

One key unknown factor that decides the neutrino flux from 
GRB host galaxies is the pion production efficiency, i.e. the 
energy loss of CR protons into pions due to collisions with 
the ISM of the host galaxy. CRs interact with the ISM along 
their path before escape, so this efficiency relates with the 
gas column density that CRs traversed. As CRs are trans¬ 
ported outward by galactic winds, the gas surface density Eg 
represents an averaged column density of matter that CRs 
traversed. However, for most GRBs, we do not have mea¬ 
surements of Eg at the GRB explosion site as GRB hosts are 
hardly spatially resolved. We note that the X-ray absorption 
column density Ah, inferred from the X-ray observations of 
the prompt and afterglow emission, has a similar role, but Ah 
reflects the column density along the direction of the line of 
sight (rather than the averaged gas column density). More¬ 
over, the value of Ah is derived assuming solar metallicity for 
the absorbing gas, while GRBs seem to occur preferentially 
in low metallicity galaxies (Stanek et al. 2006; Prieto et al. 
2008). The X-ray absorption, as measured with Swift/XRT 
observations, mainly takes place through inner shell electrons 
of metals, thus it is linked to the metallicity of the ISM in the 
host galaxies. If GRB host galaxies have lower metallicity, the 
true absorbing gas column density should be higher than the 
quoted values obtained from fitting X-ray afterglows, consid¬ 
ering this metallicity difference. We thus first study how well 
Ah can trace the gas surface density Eg in §2. We find that Ah 
is well correlated with the surface star-formation rate (Esfr) 
of the host galaxies. Then taking into account the Kennicutt 
- Schmidt law that relates the surface gas density and the sur¬ 
face star-formation rate (Kennicutt 1998), we obtain a relation 
between Ah and Eg. Once Eg is known for each GRB, we can 
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calculate the neutrino production efficiency and further calcu¬ 
late the accumulated neutrino flux from all GRB host galaxies 

(§3). 

2. THE RELATIONS BETWEEN Ah AND Eq 

Because we have very few measurements of Ssfr, but have 
measurement of A^h for many GRBs, we study whether A^h are 
correlated with Esfr and act as an indicator of the gas surface 
density. We collect the GRB hosts that have spatially resolved 
observations from literatures (Svensson et al. 2010; Fruchter 
et al. 2006; Savaglio et al. 2009), which are given in Table 
1. The total SFR and 80% light radius of these hosts Rgo are 
obtained in these observations. By comparing the 80% light 
radius and half light radius of some GRB hosts that are spa¬ 
tially resolved (Fruchter et al. 2006; Kelly et al. 2014), we 
find that Rgo — '^■Rso- Converting 80% light radius to the half- 
light radius R^q with Rgo = 2 R 50 , we calculate the surface star- 
formation rates Ssfr- With a sample of 18 GRBs that have 
both measurements of A^h from the XRT data of Swift ob¬ 
servations and Esfr of their host galaxies, we then study the 
relation between A^h and Esfr. Through the ordinary least- 
squares (OLS) bisector fitting' (Isobe et al. 1990), we find 
that 

/ \ 1.42±0.25 

Esfr = Mgyr-'kpc-^ (1) 


with a correlation coefficient of r = 0.756 and a null hypoth¬ 
esis probability of p = 7.09 x lO"”*, indicating a good correla¬ 
tion. The result of the fit is shown in the left panel of Fig.l. 
Considering the Kennicutt-Schmidt law between the surface 
SFR and the surface density of molecular plus atomic gas 
(Kennicutt 1998), i.e.. 


:(2.5±0.7)x 10” 


IM 0 PC 


M 0 yr 'kpc 


we get 

E„ = (2.1±1.0)x 10^ 


Nh 

10 ^'cm “2 


where mu is the mass of the hydrogen atom. The roughly lin¬ 
ear relation between Eg and A^h suggests that A^h is a good in¬ 
dicator of the gas surface density Eg. Note that A^h is derived 
assuming solar metallicity for the absorbing gas. Since GRB 
host galaxies usually have lower metallicity, the true absorb¬ 
ing gas column density should be corrected by the metallicity 
effect. The coefficient 2.1 ± 1.0 may reflect such a correction. 

We also study whether A^h are correlated with the total SFRs 
of the host galaxies. We collect 35 GRBs in total that have 
measurements of A^h and total SFRs, which are also listed in 
Table 1. Through the OLS bisector fitting of the data of A^h 
and SFR, we find that 

/ at \ l-55±0.26 

SFR = l 0-0 «±0'23 ( » ) M 0 yr-', (4) 

VIO cm“2/ 


with a correlation coefficient of r = 0.5576 and a null hypothe¬ 
sis probability p = 0.0014, which indicates a tight positive cor¬ 
relation. The result of the fitting is shown in the right panel of 
Fig. 1. With this relation, we can infer the SFR of GRB hosts 
with A^h for each GRB in our sample, which will be used to 
compute the galactic wind velocity in §3.1. 


* Since our goal is to estimate the functional relation between Ah and 
SsFR, and it is not clear which variable is causative. Therefore, the OLS 
bisector fitting which treat the variables symmetrically should be used, ac¬ 
cording to Isobe et al. (1990). 


3. NEUTRINO ELUX FROM GRB HOST GALAXIES 
3.1. The pion production efficiency 

The GRB accelerated CR protons travel through the 
host galaxy and produce high-energy neutrinos via proton- 
proton(f>p)-collisions with its interstellar medium (ISM). The 
collisions produce charged pions, which decay to neutrinos 
(tt"^ — 1 - , TT~ —!■ v^v^Vee~). Meanwhile, CRs can es¬ 

cape out of the galaxy through diffusion or galactic wind ad- 
vection. These two competing processes (i.e. collision and 
escape) regulate the efficiency of the pion-production of CRs, 
which can be described by /„• = 1 - exp (-fesc/Goss), where 
floss is the energy-loss time of CRs via {pp) collisions and 
fesc is the escape time of CRs. The pp-collision energy loss 

time is fioss = (nnappc) = [nYigappc/l) , where k ~ 0.5 
is the inelasticity, n is the gas number density, I is the scale 
height of the gaseous disc of the galaxy, app is the pp- 
collision inelastic cross-section, which slightly increases with 
the proton energy, given by app = (34.3-1- 1.88L-l-0.25L^) x 
r 41 ^ 

1 - () mb where L = log (Kelner et al. 

2006). Thus, the collision energy loss time is fioss = 5.4 x 

10®vr—1— { ) ( ^i>i> 

^ SOOpc y O.Olgcm-^ J VlOOmb/ 

The CRs escape from the GRB host galaxy in two 
ways. One is advective transportation by the galactic wind. 
Galactic-scale gaseous winds are ubiquitous in star-forming 
galaxies at all cosmic epochs (Heckman et al. 1990; Pettini 
et al. 2001; Shapley et al. 2003). Such winds can be driven 
by stellar winds, supernova explosions or other processes. It 
is found that, for luminous and ultraluminous infrared galax¬ 
ies at low redshifts, the winds from more luminous starbursts 
have higher speeds, roughly as oc SFR'^ with k being in the 
range of 0.25-0.35 (Martin 2005; Rupke et al. 2005; Arribas 
et al. 2014). A similar trend is found for star-forming galaxies 
at z ^ 1 (Weiner et al. 2009), i.e. 


Tw 


/ SFR \ 


(5) 


with the error in being 35%. Such type of relation be¬ 
tween the velocity of the outflowing material and SFR is ex¬ 
pected if the mechanical energy is supplied by stellar winds 
and supernova explosions. Thus, one can calculate v„ for 
each GRB hosts with measured Ah under the help of Eq. 4 
and Eq. 5. Taking v„ = 500 km s~' as the reference value for 
GRB host galaxies, the advective escape time is fadv = llv„ = 

9x10V50^(500^)"‘- 

The other way that CRs escape is diffusion, i.e., CRs are 
scattered off small-scale magnetic field inhomogeneities ran¬ 
domly and diffuse out of the host galaxy. The diffusion time 

is estimated to be Giff = jAD, where D = Dq{E/Eq) is the 
diffusion coefficient, and Dq and Eq = 3GeV are normaliza¬ 
tion factors. Since little is known about the diffusion coeffi¬ 
cient in GRB host galaxies, in the calculation we allow lower 
diffusion coefficients than that of our Galaxy, which is Dq = 
lO^^cm^s"'. The energy dependence of the diffusion coeffi¬ 
cient is also unknown and <5 = 0-1 depending on the spectrum 
of interstellar magnetic turbulence. We assume two cases, one 
is the commonly-used value <5 = 0.5, based on the measure¬ 
ment of the CR confinement time in our Galaxy (Engelmann 
et al. 1990; Webber et al. 2003). Another choice is 5 = 0.3, as¬ 
suming the Kolmogorov-type turbulence. Then the diffusive 
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escape time is fd.ff = 10 ^ 1 ( 50 ^) ( lo^w,-. )~' 
with Ep ~ 25(1 +z)s„. 

Combining the advective and diffusive escape timescales, 
the total escape time is = fdiff + Qv- When the difference 
between boss and fadv is large, fesc ~ min(fdiff,4dv)- Note that 
when fadv ^ fdiff, which holds for low-energy CRs and small 
values of Do, the pion production efficiency is independent of 
I and is given by 


( s, \ 


-1 

/ ^PP \ 

\0.01gcm 

V 500kms“* / 

UOOmb/ 


The dependence of /„• on CR energy has a break at the energy 
where fadv = 4iff- Below the break energy, /tt slightly increases 
with the energy due to the pp-collision inelastic cross-section, 
while above the break decreases as leading to a steeper 
neutrino spectrum. 

3.2. Calibration of CR flux from GRBs 

The Fermi/GBM detection rate of GRBs is about 250 per 
year (Paciesas et al. 2012), and the total gamma-ray (10-1000 
keV) fluence of GRBs detected in one year is 2 x 10“^ergcm^. 
The GBM field of view (FOV) is roughly Ogbm = 27r steradi- 
ans, so the total all-sky flux of gamma-rays in 10-1000 KeV is 
= 4 X 10“^ergcm“^yr“^. Then the energy production rate of 
gamma-rays in the energy range of 10-1000 keV, is estimated 
to be VTy(z = 0) = ^F-yH/c = 5 x 10"^^(^/0.5)ergMpc“^yr“', 
where H is the Hubble constant and ^ is a factor account¬ 
ing for the source density evolution with redshift (Eichler et 
al. 2010). According to the estimate of Katz et al. (2009), 
the present-day differential energy production rate of UHE- 
CRs above 30EeV is e\dnldep(j. = G) = (0.45 ±0.15)(a-1) x 
10^'*ergMpc“^yr“', where a is the power-law index of the CR 
spectrum. We assume a = 2 as suggested by Eermi acceler¬ 
ation. Define pp as the ratio between the differential energy 
production rate of UHECRs and the energy production rate of 
gamma-rays in 10-1000 keV. If GRBs are the source of UHE¬ 
CRs above 30 EeV, pp = (£pC/n/c/£p)/Wy(10- lOOOKeV) = 
9 ± 3, assuming f = 0.5 (Eichler et al. 2011). 

3.3. The neutrino flux 

We collect a sample of GRBs that have both measurements 
of Ah and gamma-ray fluence. As the energy coverage of 
Swift BAT is small, we choose GRBs that are detected by 
Eermi/GBM to get more accurate values of the gamma-ray 
fluence. There are 45 GRBs in total that have measurements 
of both Ah and gamma-ray fluence, as shown in Table 2. Eor 
each GRB, the (single-flavor) neutrino fluence is estimated to 
be 

2 1 2 1 

~ ~ ( 2 ) 

where is the differential energy flux of CR protons and 
Ay is the gamma-ray fluence of this GRB. The accumulated 
neutrino flux produced by all GRBs in the universe is 

^ X! 4^i^,i^GRB(47r)“' (8) 

i 

where i represents the /th GRB in our sample, A is the total 
number of GRBs in our sample, and Rgrb = 500yr“' is the 
all-sky rate of GRBs that would be detected by Eermi/GBM 
if the field of view of GBM is Att steradians. The result 
of the accumulated neutrino flux is shown in Eig.2. We 


find that the accumulated neutrino flux at PeV is (0.3 ± 
0.2) X 10“*GeV cm“^ yr“' sr“' per flavor for 6 = 0.3 and Do = 
10^®cm^s“'. The 1 sigma error in the accumulated flux is ob¬ 
tained by using formulas for the propagation of error. It results 
mainly from the uncertainties in the neutrino production effi¬ 
ciency Att, in the values of Ah, and in the CR-to-gamma-ray 
ratio rjp. 

Since the CR diffusion coefficient is not well-understood, 
we study its effect on the neutrino flux. The results for dif¬ 
ferent values of 6 and Do are shown in Pig. 3. One can see 
that, when Do decreases, the neutrino flux increases at PeV 
energies. Eor 6 = 0.3, if Do is as low as lO^^cm^s”', the neu¬ 
trino spectrum becomes too hard, although the flux increases. 
Eor 6 = 0.5, lower values of Do is needed to produce the same 
flux compared to the case of 5 = 0.3. However, the confine¬ 
ment of 100 PeV CRs producing PeV neutrinos requires that 
Do be larger than 10^^cm^s“' (Murase et al. 2013). Thus, 
even under favorable conditions about the CR diffusion coef¬ 
ficient, the neutrino flux produced by GRB CRs alone can not 
explain the IceCube observations if other assumptions used in 
our calculation are correct. 

In the calculation of the accumulated neutrino flux, we only 
considered those GRBs that triggered Eermi/GBM detector 
(i.e. adopting Rgrb = 500yr“'). There are dim GRBs that do 
not trigger the detector and their total gamma-ray fluence may 
even be larger than that of the triggered ones according to the 
simulation results in Liu & Wang (2013). These untriggered 
GRBs may also produce diffuse neutrinos as triggered GRBs. 
However, one should note that considering the contribution 
by untriggered GRBs would not change our result about the 
accumulated neutrino flux, because the the CR-to-gamma-ray 
ratio rjp needs to be re-calibrated correspondingly and thus the 
accumulated neutrino flux remains unchanged. 

There could be an exceptional case, i.e., if some dim GRBs 
do not accelerate protons to energy 10'® eV (although nu¬ 
clei can still be accelerated to ultra high energies), but they 
can still contribute to PeV neutrinos with CRs of > lOOPeV. 
Note also that this does not exclude the possibility that low- 
luminosity GRBs themselves are responsible for their origin 
(e.g., Murase & loka 2013, Bhattacharya et al. 2014). 

4. DISCUSSIONS AND CONCLUSIONS 

The above calculation has some uncertainties in the follow¬ 
ing aspects. Pirst, the estimate of the energy production rate 
of gamma-rays VEy has uncertainty (e.g., Dermer 2012). Note 
that the factor f = 0.5 is obtained by assuming that GRB rate 
follows the cosmic SPR of Porciani & Madau (2001) or that 
of Hopkins & Beacom (2008). If GRB density evolves faster 
than the cosmic SPR, f is smaller and then r/p is larger. Sec¬ 
ond, the estimate of the UHECR energy budget has uncer¬ 
tainty. The value can be a bit larger if one uses the Telescope 
Array or HiRes data. Third, the power-law index a could be 
softer than 2, then the energy budget of CRs at 100 PeV could 
be larger. Thus, given these uncertainties, the neutrino flux in 
the optimistic case could reach the observed value of IceCube. 

On the other hand, if one consider that GRB host galax¬ 
ies that are not detected by optical observations are possibly 
smaller galaxies, such as dwarf galaxies, the pion produc¬ 
tion efficiency may be smaller and thus the neutrino flux con¬ 
tributed by these galaxies would be smaller. Also, as shown 
by Eq. 6, the pion production efficiency depends on the speeds 
of the galactic winds. The properties of the galactic winds 
in GRB host galaxies are not well-explored. X-ray observa¬ 
tions of starburst galactic winds, such as M82, usually give a 



4 


higher speeds than that inferred from the optical observations 
(Strickland & Heckman 2009). If the galactic wind speeds 
of GRB hosts are proven to be higher, the neutrino flux pro¬ 
duced by GRB hosts would decrease. The pion production 
efficiency also depends on the diffusive coefficient. If Do is 
larger than lO^^cm^s"' (for example. Do = lO^^cm^s”' for our 
Galaxy), the accumulated neutrino flex would be lower than 
10“®GeV cm“^ yr“' sr“' at IPeV, as shown in Fig. 3. 

In summary, we calculated the neutrino flux produced by 
CRs accelerated by GRBs while they are propagating in the 
host galaxies based on our current knowledge about GRB and 
their host galaxies. These CRs collide with nuclei of ISM and 
produce neutrinos before they escape out of the galaxy. When 
the flux of CRs produced by GRBs is normalized with the 
observed flux of UHECRs above ^ lO'^eV, the accumulated 


neutrino flux is (0.3 ±0.2) x 10“*^GeV cm“^ s“' sr“' per flavor 
under the usual assumptions about the GRB properties and 
favorable assumptions about the CR diffusion coefficient. The 
estimate, however, has uncertainty due to uncertainty in our 
current knowledge of GRB and their host galaxies and the 
accumulated neutrino flux could reach the observed value by 
IceCube in the optimistic case, so we can not rule out this 
scenario at present. 
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Fig. 1.— Left panel: the relation between the surface star-formation rate Ssfr and the X-ray absorption column density A^h- Right panel: the relation between 
the total star-formation rate SFR and the X-ray absorption column density A^h- The data points are from table 1. The lines represent the best OLS bisector fits. 
See the text for more details. 
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Fig. 2.— The energy spectra of neutrino emission in comparison with IceCube data. The solid line represents the neutrino flux obtained using <5 = 0.3 and 
Do = lO^^cm^ s“^. and the dashed lines are the 1 -cr error. The IceCube data are also shown. 
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Fig. 3.— The energy spectra of neutrino emission for different values of S and Dq. 
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TABLE 1 

Galaxy parameters of GRB hosts 


GRB 

SFR [Mq yr-‘] 

Nh [10^‘cm--] 

Rgo [kpc] 

Ref. 

GRB 050826 

9.13 


- 

[1,2] 

GRB 051006 

51 

6+1 

- 

[3,A] 

GRB 060814 

256 

20.9+11 

- 

[3,4] 

GRB 061121 

160 

5 4 + 0.8 

- 

[3,4] 

GRB 061126 

2.38 

1 5-^-2 

- 

[1,A] 

GRB 070306 

143 

26.S^l 

- 

[3,4] 

GRB 080325 

12.9 


- 

[5,A] 

GRB 080607 

19.1 

22.8+4-’ 

- 

[5,4] 

GRB 081109 

49.0 

11 + 1.3 

^^- 1.2 

- 

[5,6] 

GRB 090328 

4.8 

0 9+1-^ 

^•^- 0.9 

- 

[7,A] 

GRB 090424 

0.8 

4 1 + 0.6 
^•^- 0.5 

- 

[8,4] 

GRB 091127 

0.22 

0.76+“-^ 

- 

[9,4] 

GRB 120624B 

23 

97+46 

26 

- 

[10,11] 

GRB 130427 

0.9 

1 

^•^- 0.1 

- 

[10,12] 

GRB 970228 

0.53 

3^ 

3.2 

[1,13,14] 

GRB 970508 

1.14 

6+'° 

1.48 

[1,13,14] 

GRB 970828 

0.87 

7.08+3:21 

2.8 

[1,15,14] 

GRB 971214 

11.4 

200 ! 2 ’« 

2.36 

[1,13,14] 

GRB 980703 

16.57 


2.42 

[1,13,14] 

GRB 990123 

5.72 


5.01 

[1,13,14] 

GRB 990705 

6.96 

13 

9.38 

[1,16,14] 

GRB 000926 

2.28 


10.25 

[1,17,14] 

GRB 011211 

4.90 

3.o!l:l 

2.69 

[1,18,14] 

GRB 020405 

3.74 

A 7 + 3.7 
'- 3.7 

11.96 

[1,19,14] 

GRB 030323 

1 

7 9 + 1-4 
'•^- 1.2 

1.86 

[20,20,14] 

GRB 041006 

0.34 

o 9 + 0.16 
^•^- 0.16 

5.19 

[1,17,14] 

GRB 050416A 

1.77 

^ 4 - K ).9 

^•^- 0.5 

2.12 

[21,4,21] 

GRB 050525A 

0.15 


1.76 

[21,4,21] 

GRB 051022 

36.46 

83.9+8;i 

2.23 

[1,5,21] 

GRB 060218 

0.05 

5 2+°-3 

_4-4-n..s_ 

0.55 

[1,2,21] 


(l)Savaglio, S., et al. 2009, ApJ, 691, 182 (2)Campana, S., et al. 2010, MNRAS, 402, 2429 (3)Perley, D. A., et al. 2014, arXiv: 1407.4456 (4)Campana, S., et al. 
2012, MNRAS, 421, 1697 (5)Perley, D. A., et al. 2013, ApJ, 778, 128 (6)Kruhler, T., et al. 2011, A&A, 534, AA108 (7)McBreen, S., et al. 2010, A&A, 516, 
AA71 (8)Jin, Z.-P, et al. 2013, ApJ, 774, 114 (9)Vergam, S. D., et al. 2011, A&A, 535, AA127 (lO)Wang, F. Y., & Dai, Z. G. 2014, ApJS, 213, 15 (ll)de 
Ugarte Postigo, A., et al. 2013, A&A, 557, LL18 (12)Littlejohns, O. M., et al. 2014, arXiv:1412.6530 (13)Galama, T. J., & Wijers, R. A. M. J. 2001, ApJ, 549, 
L209 (14)Svensson, K. M., 2010, MNRAS, 405, 57 (15) Yoshida, A., et al. 2001, ApJ, 557, L27 (16)Amati, L., et al. 2000, Science, 290, 953 (17)Kann, D. A., 
et al. 2006, ApJ, 641, 993 (18)Reeves, J. N., et al. 2003, A&A, 403, 463 (19)Covino, S., et al. 2003, A&A, 400, L9 (20)Vreeswijk, P. M., et al. 2004, A&A, 419, 

927 (21)Frachter, A. S., et al. 2006, Nature, 441, 463 (A) is from the GCN website. 
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TABLE 2 

The X-ray absorption column density and 10-1 OOOKeV fluence 


GRB 

Nh “ [ 10^' cm-^] 

GBM Fluence*^ [erg/cm^] 


141004A 

1+|.S 

1.18 X lO^'’ 

0.57 

140907A 


6.45 X 10^® 

1.21 

140703A 

5-41,4 

7.62 X 10^® 

3.14 

140512A 

0 - 521:1 

2.93 X 10-3 

0.725 

140506A 

1 - 81 :^ 

6.59 X 10-3 

0.889 

140423A 

6.6+P 

1.81 X 10-3 

3.26 

140304A 

0.381:^4 

2.24 X 10-3 

5.28 

140213A 

1 - 41.9 

2.12 X 10-3 

1.2076 

140206A 

10 - 11:3 

1.55 X 10-3 

2.74 

131105A 

i. / ^ 2 

2.38 X 10-3 

1.686 

130612A 

7.71-® 

6.80 X 10-3 

2.006 

130610A 

2 - 61 :^ 

3.54 X 10-3 

2.092 

130420A 

0.581:1 

1.16 X 10-3 

1.297 

121211A 

4 -CQ+l.O? 

-0.99 

6.41 X 10-3 

1.023 

121128A 

111 

9.30 X 10-3 

2.2 

120922A 

0.236l;P 

8.21 X 10-3 

3.1 

120907A 

-0.29 

8.09 X 10-3 

0.97 

120811C 

0 67+0-34 
-0.31 

3.45 X 10-3 

2.67 

120729A 

^•'-0.4 

5.08 X 10-3 

0.8 

120712A 

Q 14+0.31 
^•^^-0.14 

4.43 X 10-3 

4.1 

120326A 

Q 00+0.23 
^•■^^-0.22 

3.26 X 10-3 

1.798 

120119A 

1 O+0.5 
^•■^-0.4 

3.87 X 10-3 

1.73 

111228A 

9 c 1 +0.49 
^•-^^-0.47 

1.81 X 10-3 

0.715 

111107A 

0 c+6 

•^•^- 3.5 

9.07 X 10-3 

2.893 

110818A 

-0.24 

5.15 X 10-3 

3.36 

110213A 

-0.23 

9.37 X 10-3 

1.46 

110128A 

~ 0 

1.43 X 10-3 

2.339 

101219B 

0 - 51:1 

3.99 X 10-3 

0.5519 

100906A 

8 - 11:5 

2.33 X 10-3 

1.727 

100816A 

2.61:1 

3.65 X 10-3 

0.8039 

100814A 

1.6 

1.49 X 10-3 

1.44 

100728B 

2.5 

3.34 X 10-3 

2.1 

100728A 

^■'- 0.3 

1.28 X lO-"* 

1.567 

100704A 

2 1+0.9 
^•^-0.8 

1.04 X 10-3 

3.6 

100615A 

11+1.3 

8.72 X 10-3 

1.398 

091208B 

8.31-3 

6.19 X 10-3 

1.063 

091127 

0.761-1 

2.07 X 10-3 

0.49 

091020 

5 - 81:6 

8.35 X 10-3 

1.71 

090926B 

13.9+1-S 

1.08 X 10-3 

1.24 

090516A 

22.91, 

1.72 X 10-3 

4.1 

090424 

4 1+0.6 
^•^-0.5 

4.63 X 10-3 

0.544 

090423 

102l« 

8.16 X 10-'' 

8 

081221 

26-11:1 

3.00 X 10-3 

2.26 

080916A 

0+3.2 
° 1.9 

7.81 X 10-3 

0.689 

080905B 

22.61-3 

2.91 X 10-3 

2.374 


(a)From the GCN website. (b)From Fermi-GBM GRB list of detections, http://www.asdc.asi.it/grbgbm/ 





